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We report the first measurement of transverse single-spin asymmetries in J/tp production from 
transversely polarized p + p collisions at y/s = 200 GeV with data taken by the PHENIX experiment 
in 2006 and 2008. The measurement was performed over the rapidity ranges 1.2 < \y\ < 2.2 and \y\ < 
0.35 for transverse momenta up to 6 GeV/c. J/tp production at RHIC is dominated by processes 
involving initial-state gluons, and transverse single-spin asymmetries of the J/tp can provide access 
to gluon dynamics within the nucleon. Such asymmetries may also shed light on the long-standing 
question in QCD of the J/tp production mechanism. Asymmetries were obtained as a function of 
J/tp transverse momentum and Feynman-x, with a value of -0.086 ± 0.026 stat ± 0.003 syst in the 
forward region. This result suggests possible nonzero trigluon correlation functions in transversely 
polarized protons and, if well defined in this reaction, a nonzero gluon Sivers distribution function. 

PACS numbers: 14.20.Dh, 13.88.+C, 13.85.Ni, 14.40.Pq 



I. INTRODUCTION 

The transverse single-spin asymmetry (SSA) quanti- 
fies the asymmetry of particle production relative to the 
plane defined by the transverse spin axis and the momen- 
tum direction of a polarized hadron. SSAs have come to 
be recognized as a means of accessing QCD dynamics, 
both within initial-state hadrons and in the process of 
hadronization from partons. Large azimuthal transverse 
single-spin asymmetries of up to ~ 40% were first ob- 
served at large Feynman-x (xp = 2-ppj^fs, where pl is 
the momentum along the beam direction) in pion pro- 
duction from transversely polarized p + p collisions at 
\/s = 4.9 GeV in 1976 l], contrary to theoretical ex- 
pectations at the time Q, and subsequently observed 
in hadronic collisions over a range of energies extend- 
ing up to y/s — 200 GeV d-Q. In recent years nu- 
merous measurements of transverse SSAs have been per- 
formed in semi-inclusive deep-inelastic scattering (SIDIS) 
off a transversely polarized proton or deuteron target 
as well jT^-[l^ . In order to describe the large trans- 
verse SSAs observed, two approaches have been devel- 
oped since the 1990s, after early pioneering work by Efrc- 
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mov and Teryaev [161 Il7j . One approach requires higher- 
twist contributions in the collinear factorization scheme, 
i.e. the exchange of a gluon between one of the partons 
taking part in the hard scattering and the color field of 
either an initial- or final-state hadron. This was first pro- 
posed by Qiu and Sterman for gluon exchange in the ini- 
tial state 18[ and by Kanazawa and Koike for exchange 
in the final state Gluon exchange in either the ini- 
tial or final state leads to terms including multiparton 
correlation functions, which can generate a nonzero SSA. 
The other approach utilizes parton distribution functions 
and/or fragmentation functions that are unintegrated in 
the partonic transverse momentum, kx; these functions 
are generally known as transverse-momentum-dependent 
distributions (TMDs). These two approaches have dif- 
ferent but overlapping kinematic regimes of applicability 
and have been shown to correspond exactly in their re- 
gion of overlap [20[ ■ 

While higher-twist parton-nucleon spin-momentum 
correlation functions and TMDs were born within the 
nucleon structure community, as understanding of them 
matures, their implications for other areas of QCD are 
starting to be realized. J/ip production has been ex- 
tensively studied over the last decades, but the details 
of the production mechanism remain an open question 
(see ref. [2lj for a comprehensive review), and an addi- 
tional complication is that approximately 30%-40% of 
the measured J/tp mesons in hadronic collisions are pro- 
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duced indirectly from feed-down from ip' and Xc |22| . It 
was proposed in 2008 by Yuan [23| that within the frame- 
work of nonrelativistic QCD (NRQCD) [24j, the trans- 
verse SSA of J/ip production can be sensitive to the J/ip 
production mechanism, assuming a nonzero gluon Sivers 
function [Hj], which is a TMD that describes the corre- 
lation between the transverse spin of the proton and the 
kx of the partons within it. Specifically, Yuan predicts 
that a nonzero gluon Sivers function will produce a fi- 
nite transverse SSA for color-singlet J / ip production (2(| 
in p + p collisions, but the asymmetry should vanish for 
color-octet production 24] in p + p due to cancelation 
between initial- and final-state effects, while a nonzero 
asymmetry for J / ip production in SIDIS is only expected 
within the color-octet model. It should be noted that 
the relationship between the transverse SSA and the pro- 
duction mechanism is not quite as simple in the collinear 
higher-twist approach, with partial but not full cance- 
lation of terms [13] in the cases where the asymme- 
try uniformly vanishes in the TMD approach presented 
by Yuan. Another important point to note regarding 
the TMD as compared to the collinear, higher-twist ap- 
proach is that very recent theoretical work [28( suggests 
that factorization of hard processes in perturbative QCD 
(pQCD) into transverse-momentum-dependent distribu- 
tion and fragmentation functions convoluted with par- 
tonic hard-scattering cross sections is not valid for pro- 
cesses involving more than two hadrons. Thus, in the 
process p+p — > J/ip + X a gluon Sivers function may not 
be well defined; however, the definition within a factor- 
ized pQCD framework of the corresponding trigluon cor- 
relation functions in the collinear, higher-twist approach 
is believed to be valid. 

Measurements of heavy flavor transverse SSAs in p + p 
collisions are of interest because they serve to isolate 
gluon dynamics within the nucleon, with heavy quarks 
being dominantly produced via gluon-gluon interactions. 
Very little is presently known about trigluon correlation 
functions or gluon TMDs. Measurements of the trans- 
verse momentum (pt) spectrum of bottomonium produc- 
tion have been used to extract the fcy-unintegrated distri- 
bution of gluons in an unpolarized proton [29j . Similarly, 
measurements of open heavy flavor transverse SSAs have 
been proposed as a means to isolate gluon TMDs and/or 
their corresponding twist-thr ee g luon correlation func- 
tions in polarized protons |30l - [32l [. A previous PHENIX 
measurement of the transverse SSA in neutral pion pro- 
duction at midrapidity [33j j as well as measurements by 
the COMPASS collaboration of the SSA in semi-inclusive 
charged hadron production [Til ], both consistent with 
zero, have been used to provide initial constraints on the 
gluon Sivers function and suggest that it is small [13, [Hj] , 
but further data relevant to gluon TMDs are greatly 
needed. 

In this paper the first measurement of transverse SSAs 
in J/ip production is presented. The data were taken 
by the PHENIX experiment at the Relativistic Heavy 
Ion Collider (RHIC) during the 2006 and 2008 polarized 



proton runs at y/s — 200 GeV. The pr and xp dependen- 
cies are studied, for rapidity regions of— 2.2 < y < —1.2, 
\y\ < 0.35, and 1.2 < y < 2.2, and p T up to 6 GeV/c. 



II. ANALYSIS 



A. Measuring Transverse Single-Spin Asymmetries 



Transverse single-spin asymmetries lead to modula- 
tions of the cross section in the azimuthal angle due to 
the projection of the polarization vector into the direc- 
tion of the produced particle. Our measurement is per- 
formed in two separate hemispheres referred to as 'left' 
and 'right' where left is defined as the axis which forms a 
right-handed coordinate system with the beam momen- 
tum vector and one of the spin orientations, denoted as 
f. For a vector S in the direction of the f spin and beam 
momentum P, left is defined as p ■ (^S x pj > with p 

being the momentum vector of the outgoing particle. 

The left-right transverse SSA can be extracted using 
Eq. ([1} . This equation applies to particle yields observed 
to the left side of the polarized beam. 



A 



N 



V (at + <A) ' 



(1) 



where a^(a^) represents the production cross section 
with beam polarized in the f (4) direction, integrated 
over the left hemisphere, and V is the beam polarization. 
An overall minus sign is required for An on the right side 
of the polarized beam. 

The geometric scale factor / corrects for the convo- 
lution of an azimuthal asymmetry with detector accep- 
tance. For a sinusoidal asymmetry, as generated by the 
Sivers function, the factor becomes 



/ = 



f£ e(<p) sin 4>d(j) 
Jq e{(p)d<p 



(2) 



where <p is the azimuthal angle between the outgoing par- 
ticle and the proton spin, and e{4>) is the efficiency for 
detecting a J/ip at a given <f>. The limits of integration 
correspond to the hemisphere in which the measurement 
is being made. It should be noted that even for a detector 
with full azimuthal coverage the factor / is not unity be- 
cause we are not determining the amplitude of the mod- 
ulation directly but measuring an asymmetry which is 
integrated over entire hemispheres in azimuth. 

Both proton beams in RHIC were polarized. In order 
to derive single-spin asymmetries, we sum over polariza- 
tion direction states in one of the two beams: 



(3) 
(4) 
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where a again indicates the ^-dependent cross section in- 
tegrated over one hemisphere, and f and 4- represent the 
spin orientations of the two beams. Differences in lumi- 
nosities can lead to false asymmetries; see Section III Dl 
We perform the analysis separately for each beam. Up 
to 111 out of a possible 120 bunches were filled in each 
RHIC ring, with 106 ns between bunches. Four pre- 
set spin patterns with approximately equal numbers of 
bunches polarized in opposite directions were alternated 
during the transverse running periods to minimize sys- 
tematic effects from the injection process. Stores were 
nominally held for eight hours. The spin orientation at 
the PHENIX interaction point was maintained at the de- 
fault vertical orientation in the RHIC ring for 2008 data 
taking and rotated to radial orientation for 2006. 

Explicitly including the spin orientations of both 
beams, Eq. ((T|) can be rewritten as: 



.4 



JY 



V (crtt + a U) + ( CT 4.t + a ii) 



V (ATtt + TZxNU) + (K 2 N^ + K 3 NU) 



(5) 



where N^, N^, and are the experimental yields 
in each spin configuration, and 1Z± = / C^^, 1Z 2 = 
£tt J £-lt a nd 7^ 3 — £tt J £ll are ra tios of the pro- 
vided luminosities C in each spin orientation. The use 
of a single polarization factor assumes that bunch po- 
larizations V are equal regardless of spin orientation in 
each beam. We have also assumed that detector efficien- 
cies are equal for each spin orientation, given that the 
orientations change every 106 or 212 ns. 

Beam polarizations at RHIC are measured by two dif- 
ferent polarimeters, a fast carbon target polarimeter [36[ 
for relative polarization measurements and a hydrogen 
jet polarimeter [37j for an absolute measurement. The 
carbon target polarimeter is used to measure relative 
beam polarization store-by-store, and the hydrogen jet 
polarimeter is used for absolute calibration of the beam 
polarization. During the 2006 run the average transverse 
beam polarizations were 

0.53 ± 0.02 s y st (clockwise) 

0.52 ± 0.02 s y st (counterclockwise). 

Fill-to-fill variation of beam polarization was ±0.03 
(1 a) for the clockwise beam and ±0.04 (1 a) for the 
counterclockwise beam. The labels represent the direc- 
tion in which the beam is circulating when looking from 
above, and the systematic uncertainties are uncorrelated 
between beams. There is an additional systematic uncer- 
tainty of 3.4% correlated between the two beams. The 
average beam polarizations during the 2008 data taking 
were 

0.48 ± 0.02 s y st (clockwise) 

0.41 ± 0.02 s y st (counterclockwise) 



each with a fill-to-fill variation of ±0.04 (1 a), and with 
an additional systematic uncertainty of 3.0% correlated 
between the beams. 

For statistically limited measurements it may be im- 
possible to measure asymmetries using separate yields 
based on the spin orientation of both beams. In this case 
a good approximation can be made using 



A N = 



(6) 



where the spin orientation of the polarized beam is de- 
noted by the arrow, and we have a single relative luminos- 
ity 1Z = jr£. Systematic effects which will be discussed in 
Section III D I can be reduced by bringing 1Z to a constant 
value. 

To further simplify Eq. (|6|) we can eliminate all ex- 
plicit uses of the relative luminosity and calculate a single 
asymmetry for both the left and right hemispheres using 



A 



N 



(7) 



The geometric scale factor 
e{4>) sin ^d(f) 



f' = 2 



{•7T 

Jo e 



Ce(4>)d4> 



(8) 



is different from the scale factor in Eq. and Eq. © 
because left and right are treated simultaneously, lead- 
ing to differences between the two expressions of order 



B. J/ip Measurements 

Measurements were carried out by the PHENIX ex- 
periment at RHIC, where the cross section and polariza- 
tion of J/tp mesons in y/s = 200 GeV p + p collisions 
have also been measured [HI, [3!|. An overview of the 
PHENIX apparatus can be found in [4(|. At forward 
and backward rapidities, J/tp /x 4 \x~ were measured 
with two muon spectrometers [4l|, for 1.2 < \y\ < 2.2 
and A(f> — 2ir, using data recorded in 2006 and 2008. At 



midrapidity, asymmetries were studied via J/tp — > e 
with the central arm spectrometers [42j ]. for \y\ < 0.35 
and A(p = 2 X | , using 2006 data. Collisions were identi- 
fied by triggering on a valid collision vertex as measured 
by the beam-beam counters (BBC) [H[. 

Muon track candidates were detected at forward and 
backward rapidities with respect to the polarized beam 
using two muon spectrometers on the north and south 
sides of the experiment. Each spectrometer consists of 
3 cathode strip tracking chambers in a magnetic field 
(MuTr) and 5 layers of Iarocci tube planes interleaved 
with thick steel absorbers (MuID). 

Muon candidates were selected using the BBC trigger 
in coincidence with combinations of hits in the MuID 
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called symsets. Each symset is defined by projecting from 
the center of the interaction region through a tube in the 
first layer of the MuID to a window of ±3 additional 
tubes in each subsequent layer. Muon candidates are 
characterized by the depth of their penetration through 
the MuID. A deep track requires at least one hit of the 
symset to occur in either the first or second layer and one 
hit in either the fourth or fifth layer with at least three 
layers containing hits, while a shallow track requires only 
one hit in either the first or second layer with at least two 
hits in the first three layers. In order for a muon to reach 
the MuID it must have a minimum total momentum of 
~1.9 GeV/c, and ^2.4 GeV/c is needed to reach the last 
layer. 

The dimuon trigger for the 2006 data required at least 
one muon candidate satisfying the deep trigger require- 
ment and another satisfying at least the shallow trigger, 
while the trigger for the 2008 data required two deep 
muon candidates. After run selection, the muon spec- 
trometer on the north (south) side of PHENIX sampled 
an integrated luminosity of 1.75 (1.63) pb _1 in 2006 and 
4.33 (4.30) pb" 1 in 2008. 

For candidate muon pairs, the collision vertex was re- 
quired to be within 35 cm of the center of the interaction 
region along the beam direction, and each track was re- 
quired to have longitudinal momentum 1.4 < p z (GeV/c) 
< 20. The distance between the track projections from 
the MuID and the MuTr to the first MuID layer was re- 
quired to be less than 25 (30) cm for the detector on 
the north (south) side, and the angular difference was 
required to be less than 10 degrees. A fit to the common 
vertex of the two tracks was performed and required to 
have a x 2 < 20 for 4 degrees of freedom. 

For an electron candidate at midrapidity, a coincidence 
was required between the BBC trigger and a trigger de- 
signed to select electrons (ERT), which required a nom- 
inal minimum energy of 0.4 GeV in a 2 x 2-tower re- 
gion (Ar) x A<j) = 0.02 x 0.02) of the electromagnetic 
calorimeter (EMCal) and corresponding activity in the 
ring-imaging Cerenkov detector (RICH) . After run selec- 
tion the central arm spectrometers sampled an integrated 
luminosity of 1.36 pb -1 in 2006. 

The momentum of electron candidate tracks was re- 
constructed using the drift chamber (DCH), and their 
energy was measured with the EMCal. A minimum mo- 
mentum cut of 0.5 GeV/c was applied, and the ratio of 
the measured energy to momentum (E/p) was required 
to be within four standard deviations (a) of 1, with the 
resolution in E/p varying with p between 10% and 15%. 
Position matching between the track in the DCH and 
energy cluster in the EMCal was required to be < 4a in 
both the beam direction and azimuth, and only collisions 
occurring within 30 cm of the center of the interaction re- 
gion along the beam direction were considered, matching 
the central arm acceptance. 



TABLE I: Background asymmetries as a function of px for 
PHENIX muon spectrometers. The uncertainties given are 
statistical. 



p T (GeV/c) 


< Xf > 


data set 


aBG 
A N 


0-6 


-0.081 


2006 


-0.003T0.028 




-0.082 


2008 


-0.072±0.034 




0.084 


2006 


-0.008±0.028 




0.086 


2008 


-0.003T0.035 


0-1.4 


-0.081 


2006 


-0.021±0.034 




-0.081 


2008 


-0.089±0.042 




0.085 


2006 


0.002±0.034 




0.087 


2008 


-0.008T0.043 


1.4-6 


-0.081 


2006 


0.001T0.053 




-0.082 


2008 


-0.041±0.068 




0.084 


2006 


-0.039±0.053 




0.086 


2008 


0.024±0.066 



C. Analysis Method for J ftp — > fj, + jj, 

Transverse single-spin asymmetries for the J/tp — > 
decay channel were determined by subtracting a 
background asymmetry from the inclusive signal as 



.4 



N 



\ Incl 
N 



r-A N G 



(9) 



where the Am values on the right-hand-side were cal- 
culated using Eq. ([5]). The asymmetry A 1 ^ 1 is for 
oppositely-charged muon pairs in the invariant mass 
range ±2<r around the J/ip mass (where a is the mass 
resolution of the detector), and A N G for the analysis of 
the 2006 data set is the asymmetry for oppositely-charged 
muon pairs in the invariant mass range 1.8 < m (GeV/c 2 ) 
< 2.5 along with charged pairs of the same sign in in- 
variant mass range 1.8 < m (GeV/c 2 ) < 3.6 (shaded 
areas in Fig.QJa) and (b) respectively). For the analysis 
of the 2008 data set the lower limit of the mass range 
was 2.0 GeV/c 2 . As a cross check, A N G was calculated 
separately for the oppositely-charged muon pairs and 
charged pairs of the same sign, and the results were con- 
sistent. Table U gives the measured values of A N G . The 
J/ip mass resolution was measured to be 0.143 ± 0.003 
(0.154 ± 0.004) GeV/c 2 in the muon spectrometer on the 
north (south) side for the 2006 data set and consistent 
with this for the 2008 data set. 

The background fraction r is defined as 



Nl nc l — Nj/^f, 
Nlncl 



(10) 



where Nj nc i is the total number of oppositely-charged 
muon pairs in the invariant mass range ±2<r around 
the J/ip mass, and Njm is the J/ip yield in the same 
mass range. Nj/^ was extracted by fitting the invari- 
ant mass spectrum of oppositely-charged muon pairs with 
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FIG. 1: (color online). Invariant mass spectra for (a) 
oppositely-charged muon candidate pairs and (b) charged 
candidate pairs with the same sign for the muon spectrom- 
eter on the north side of PHENIX from the 2006 data set. 
The solid line is the sum of the J/ip (dashed) and ip' (dot- 
ted) Gaussians, along with a third-order polynomial (dotted- 
dashed) background. 



two Gaussians (for the J/ip and ip' resonances, the peak 
shape of which is dominated by the detector resolution) 
and a third-order polynomial (for the remaining pairs), 
using the expression: 



(a +a\M 



N 



a 2 M 2 + a 3 M 3 ) 



[M—M. i ) 



2vr, 



2ny 



(11) 



The free parameters in the fit are the four polynomial 
parameters as well as Nji^, Mj/^, a, and N^>. The mass 
and width of the ip' are fixed relative to the J/ip based 
on simulations. The polynomial is used to fit the back- 
ground from both physical sources (i.e. Drell-Yan, open 
heavy flavor) and uncorrelated track combinations. Un- 
correlated track combinations comprise more than 50% of 
the oppositely-charged muon pairs under the J ftp mass 
peak. An example mass spectrum and fit is shown in 
Fig. [Ha). The total J ftp yield in the dimuon decay chan- 
nel was 6403 ± 126 for the 2006 and 15380 ± 150 for the 
2008 data set. Table [TT] gives background fractions in 
different px ranges determined from the fit. The asym- 
metries were measured store- by-store, using Eq. ([3]), with 
the final results obtained by averaging over all stores. 



TABLE II: Total background fractions as a function of pr for 
muon spectrometers on the north and south sides of PHENIX. 
Backgrounds were higher in the 2006 data set because the 
less restrictive trigger requirement allowed more random track 
combinations. 



p T (GeV/c) 


data set 


detector 


background fraction (%) 


0-6 


2006 
2006 
2008 


South 
North 
South 
North 


21.7±0.6 
19.1±0.4 
16.4±0.2 
14.2±0.2 


0-1.4 


2006 
2006 
2008 
2008 


South 
North 
South 
North 


23.2±0.7 
22.0±0.7 
16.1±0.3 
15.5±0.3 


1.4-6 


2006 
2006 
2008 
2008 


South 
North 
South 
North 


20.1±0.8 
14.1±0.5 
15.6±0.4 
10.5±0.2 



A geometric scale factor from 2006 data of / = 1.57 ± 
0.04 was determined from J jip azimuthal distributions in 
data and was found to be independent of px within sta- 
tistical uncertainties. In the 2008 data, the polarization 
of the clockwise circulating beam was found to be rotated 
by 0.25 ± 0.033 rad away from vertical (see Appendix p), 
meaning that / was different depending on which beam 
was considered polarized in the analysis. The geometric 
scale factors from that analysis were / = 1.64 ± 0.01 for 
the clockwise circulating beam and / = 1.56 ± 0.01 for 
the counter-clockwise circulating beam 

Several systematic checks were performed in the analy- 
sis. Asymmetries were determined using random spin ori- 
entations to cross check their consistency with zero. The 
means were found consistent with zero and widths con- 
sistent with statistical uncertainties in A^. The parity- 
odd asymmetries along the proton spin direction were 
also measured and found to be consistent with zero as 
expected for the strong interaction. 

Following previous PHENIX publications systematic 
uncertainties in this analysis are categorized as Type A, 
point-to-point uncorrelated, Type B, scaling all points 
in the same direction but not by the same factor, and 
Type C, scaling all points by the same factor. It should 
be noted that all scale uncertainties affect both the cen- 
tral values and the statistical uncertainties such that the 
statistical significance of the measurement from zero is 
preserved. A single Type A systematic uncertainty is in- 
cluded, which is due to the fit used to determine the back- 
ground. The fit was nominally performed using a range 
of 1.8 < m (GeV/c 2 ) < 5.0 for the 2006 and 2.0 < m 
(GeV/c 2 ) < 5.0 for the 2008 data. To estimate the un- 
certainty due to the choice of fit range, the fit was then 
performed using a range of 1.5 < m (GeV/c 2 ) < 5.5, 
and the difference of the calculated An with the nominal 
was taken as a systematic uncertainty. Other systematic 
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uncertainties are discussed in Section Hill 

D. Analysis Method for J/tj) — > e + e~ 

For the RHIC luminosities and store lengths in 2006 
there were too few J/ips detected in the PHENIX central 
arms to calculate a separate asymmetry for each store. 
Instead, asymmetries were calculated using Eq. ([7]) from 
statistics integrated across the four-week running period. 
The integration of statistics over a long time-period re- 
quires special care to be taken in order to avoid intro- 
ducing false asymmetries. 

If we measure an asymmetry in n time periods with, 
for simplicity, equivalent statistics in each measurement 
i, we can determine an averaged asymmetry in the left 
hemisphere using Eq. ([5]) as 

n f^ViNj + KiNf n 

Since each measurement is probing the same physical ob- 
servable we have the same A^ ti for all i within statistical 
uncertainties. This allows us to calculate the asymmetry 
as 

n 

An = ^ , (13) 

i=l 

which corresponds to integrating our n measurement pe- 
riods into a single measurement. For a statistically lim- 
ited measurement we do not calculate 1Zi, fi, and Vi for 
all i but assume that the measurement can be made as 

n 

An = ^ , (14) 

i=l 

where the brackets denote a luminosity-weighted average 
over the course of the measurement. In order for the same 
physical observable to be calculated by this expression 
without additional systematic uncertainty we must have 
TZi = (11), fi = (/), and Vi = (V) for all i. 

The PHENIX central arm acceptance was stable 
enough over the course of the running period so that the 
assumption is valid for /, but variations in the polariza- 
tion and relative luminosity must be taken into account. 
Variations in the polarization were found to be consistent 
with statistical and systematic errors on measurements of 
the polarization and contribute to an uncertainty in the 
overall scale of Am (scale uncertainties will be discussed 
in more detail in Section Mil) . Variations in 7Z have the 
most significant effect on the asymmetry as they can po- 
tentially contribute to false asymmetries. 



In order to stabilize relative luminosities a procedure 
was developed wherein several bunches of colliding pro- 
tons are removed from the analysis so that the relative 
luminosity within each store is brought as close to unity 
as possible given the finite number of bunches. In prin- 
ciple we could choose any constant value for the relative 
luminosity, but unity is chosen for convenience. 

First, all bunches with luminosities greater than 2 stan- 
dard deviations away from the mean bunch luminosity of 
a store are removed. A bunch is then chosen at random, 
and if removing this bunch from the analysis brings the 
relative luminosity closer to unity it is removed. Oth- 
erwise it is kept, and another bunch is chosen at ran- 
dom. The process continues until the relative luminosity 
is within 1% of unity or as close to unity as possible given 
the finite number of bunches. The corrected relative 
luminosities were distributed with an RMS of approxi- 
mately 1.5% away from unity, and the entire procedure 
removed approximately 5% of the provided luminosity 
from the data sample used in this analysis. 

Since removed bunches are chosen at random, the re- 
sult of the analysis is not unique, and a systematic uncer- 
tainty is introduced. In order to determine this system- 
atic uncertainty, the analysis was run 5000 times, and the 
resulting asymmetries and statistical uncertainties were 
histogrammed. The mean values of these two histograms 
were then taken as the central value and statistical un- 
certainty, and the RMS of the histogrammed asymmetry 
was taken as the dominant systematic uncertainty in the 
analysis. 

It is not generally true that detector efficiencies are 
independent of spin orientation in the PHENIX central 
arms when triggering with the ERT, as even and odd 
numbered crossings are triggered by separate circuits 
with slightly different gains. While the results for even 
and odd crossings are consistent within statistical uncer- 
tainties and no systematic trend is observed, to eliminate 
any possible effects from the different trigger circuits we 
measure asymmetries separately for even and odd cross- 
ings and combine the resulting asymmetries. 

The PHENIX central arms detected 539 ± 25 J/ip 
mesons in an invariant mass range 2.7 < m (GeV/c 2 ) 
< 3.4. The resolution of the J/ip mass peak was found 
to be 0.061 ± 0.002 GeV/c 2 . Charged lepton pairs with 
the same sign are due to uncorrelatcd tracks and are sub- 
tracted from the oppositely-charged pairs using the num- 
ber of like-signed pairs as 2y / N e + e +N e - e - . Both mass 
spectra can be found in Fig. [5J After subtracting the 
uncorrelated tracks from the mass spectra, there is a re- 
maining continuum background from open heavy flavor 
and Drell-Yan production as well as a small number of 
lepton pairs from rp' decays reconstructed to low mass. 
The contribution from conversion electrons is negligible. 
The background fraction from continuum pairs as defined 
in Eq. (jlOl) can be found in Table IIII1 There were not 
enough statistics in the continuum background to deter- 
mine A N G , and inclusion of such background was not 
found to significantly affect the signal asymmetry. An 
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TABLE III: Fraction of measured electron pairs coming from 
the continuum background in the central spectrometer. 



P t (GeV/c) 


r (%) 


0-6 


6.6±0.4 


0-1.4 


5.6±0.5 


1.4-6 


7.8±0.7 


overall dilution of the si{ 
A B N G =0. 


mal is included assuming that 


TABLE IV: Geometric scale factors determined from simula- 
tion for J/ip — > e + e~ . 


Pt (GeV/c) 


/' 


0-6 


1.62±0.01 


0-1.4 


1.61±0.01 


1.4-6 


1.70±0.02 



There were not enough J ftps detected in order to de- 
termine geometric scale factors for the asymmetries at 
midrapidity from data. These factors were instead de- 
termined using a GEANT [44j Monte Carlo simulation 
of single J/ip decays with a full geometric description 
of the central arm detectors including all known ineffi- 
ciencies. Geometric scale factors used in the analysis are 
given in Table HVl and mean polarizations in Section III Al 
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FIG. 2: (color online). Invariant mass spectra for oppositely- 
charged electron candidate pairs and uncorrelated track pairs 
defined as 2y/N e + e + N e - e - in the central spectrometer. 

A number of systematic checks were performed to en- 
sure the validity of the measurement. The number of 
J/ips produced divided by the provided beam luminosity 
was calculated and found to be constant across stores. 
Asymmetries were determined with randomized spin di- 
rections as described in Section III CI and found to be dis- 
tributed in normal distributions which had means consis- 
tent with zero and widths consistent with the statistical 



uncertainties. Parity-odd asymmetries along the proton 
spin direction were measured and found to be consistent 
with zero, as expected. 



III. RESULTS AND SUMMARY 

Figures|3]and|3]present the measured transverse single- 
spin asymmetry in J/ip production versus xf, and Fig. [5] 
shows the measured transverse SSA at different rapidi- 
ties as a function of pr- The results are tabulated in 
Table [V] The Type A systematic uncertainties, which 
are point-to-point uncorrelated, are due to the procedure 
used to stabilize the relative luminosity in the midrapid- 
ity analysis, as described in Section III Dl and the fit used 
to determine the background in the analysis at forward 
and backward rapidities, as described in Section III CI As 
discussed in Section Hi Al uncertainties in the geometric 
scale factors and polarizations lead to a fractional scale 
uncertainty on An- The uncertainties shown in the last 
two columns of Table IVl are Type B systematic uncertain- 
ties, scaling all xf or pt points in the same data set in 
the same direction but not necessarily by the same factor. 
There are additional Type C scale uncertainties, which 
scale all points in exactly the same way, due to the fully 
correlated polarization uncertainties in each data set of 
3.4% for 2006, 3.0% for 2008, and 2.4% for the combined 
2006 and 2008 data sets. It should be noted that all 
scale uncertainties affect both the central values and the 
statistical uncertainties such that the statistical signifi- 
cance of the measurement from zero is preserved. While 
in principle the polarization uncertainties do not affect 
An symmetrically due to the fact that An oc -p, the dif- 
ference in the value of the uncertainties scaling to larger 
and smaller magnitudes of An is less than the precision 
shown. 

As the functional form of the asymmetry in if and 
Pt is completely unknown, no correction has been made 
for potential smearing effects. A simulation study was 
performed assuming a linear dependence of on if, 
and it was found that smearing effects were less than 10% 
of the value of the input asymmetry. 

The measured asymmetry at forward xf is negative, 
-0.086 ± 0.026 ± 0.003, with a statistical significance 
from zero of 3.3<r, suggesting a nonzero trigluon cor- 
relation function in transversely polarized protons and, 
if well defined as a universal function in the reaction 
p + p — > J/tp + X, a nonzero gluon Sivers function. A 
nonzero transverse SSA in J /%p production in p + p gen- 
erated by gluon dynamics may seem surprising given the 
SSAs consistent with zero in midrapidity neutral pion 
production at PHENIX [33] and semi-inclusive charged 
hadron production at COMPASS [ll|. However, the de- 
tails of color interactions have been shown to play a ma- 
jor role in SSAs [28| . so further theoretical development 
will be necessary before we fully understand the relation- 
ships among these measured asymmetries. As discussed 
in Ref. [231 ] . a nonzero transverse SSA in J/%p produc- 
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TABLE V: An vs. pr in forward, backward and midrapidity. Systematic uncertainties in the last two columns are due to the 
geometric scale factor and the polarization, respectively. There are additional Type C uncertainties due to the polarization of 
3.4%, 3.0%, and 2.4% for the 2006, 2008, and combined 2006 and 2008 results. See text for details. 



Pt 


Data Sample 


< x F > 


A N 


8A N 


SA N 


SA f N (%) 


SA N (%) 


(GeV/c) 








(stat) 


(Type A syst) 


(Type B syst) 


(Type B syst.) 




2006 


-0.081 


-0.024 


0.044 


0.003 


0.6 


2.3 




2008 


-0.082 


-0.010 


0.032 


0.004 


0.4 


3.4 




2006 + 2008 


-0.081 


-0.015 


0.026 


0.002 


0.4 


2.8 


0-6 


2006 


0.000 


-0.064 


0.106 


0.026 


0.6 


2.3 




2006 


0.084 


-0.105 


0.044 


0.005 


0.6 


2.3 




2008 


0.086 


-0.075 


0.032 


0.003 


0.4 


3.3 




2006 + 2008 


0.085 


-0.086 


0.026 


0.003 


0.4 


2.7 




2006 


-0.081 


0.050 


0.067 


0.007 


0.6 


2.3 




2008 


-0.081 


-0.025 


0.046 


0.008 


0.4 


3.4 




2006 + 2008 


-0.081 


-0.001 


0.038 


0.005 


0.4 


2.8 


0-1.4 


2006 


0.000 


-0.063 


0.128 


0.031 


0.6 


2.3 




2006 


0.085 


-0.065 


0.066 


0.005 


0.6 


2.3 




2008 


0.087 


-0.064 


0.045 


0.003 


0.4 


3.4 




2006 + 2008 


0.086 


-0.064 


0.037 


0.003 


0.4 


2.7 




2006 


-0.081 


-0.073 


0.065 


0.002 


0.6 


2.3 




2008 


-0.082 


-0.023 


0.046 


0.010 


0.4 


3.5 




2006 + 2008 


-0.082 


-0.039 


0.038 


0.002 


0.4 


2.8 


1.4-6 


2006 


0.000 


-0.068 


0.188 


0.045 


1.2 


2.3 




2006 


0.084 


-0.046 


0.064 


0.005 


0.6 


2.3 




2008 


0.086 


-0.073 


0.046 


0.007 


0.4 


3.3 




2006 + 2008 


0.085 


-0.064 


0.037 


0.004 


0.4 


2.7 



0.2 
0.15 

0.1 
0.05 

* -° 

-0.05 
■0.1 
-0.15 



- p+p -> J/y+X a\\fs = 
: <P T >=1 .6 GeV/c (side po 
-_ <p T >=1 .5 GeV/c (middle 


200 GeV 
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FIG. 3: (color online) Transverse single-spin asymmetry in 
J/ip production as a function of xf for 2006 and 2008 data 
sets. The error bars shown are statistical and Type A system- 
atic uncertainties added in quadrature. Type B systematic 
uncertainties are not included but are 0.003 or less in abso- 
lute magnitude and can be found in Table [Vj An additional 
uncertainty in the scale of the ordinate due to correlated po- 
larization uncertainties of 3.4% (3.0%) for the 2006 (2008) 
data set is not shown. See text for details. 



FIG. 4: (color online) Transverse single-spin asymmetry in 
J/ip production as a function of xf for combined 2006 and 
2008 data sets. The error bars shown are statistical and Type 
A systematic uncertainties added in quadrature. Type B sys- 
tematic uncertainties are not included but are 0.002 or less in 
absolute magnitude and can be found in Table [V] An addi- 
tional uncertainty in the scale of the ordinate due to correlated 
polarization uncertainties of 2.4% (3.4%) for the points with 
\xf\ > (xf = 0) is not shown. See text for details. 
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FIG. 5: (color online) Transverse single-spin asymmetry of 
J/ip mesons plotted against J/ip transverse momentum. See 
Table IVI for mean xf values for each point. The error bars 
shown are statistical and Type A systematic uncertainties 
added in quadrature. Type B systematic uncertainties are 
not included but are 0.002 or less in absolute magnitude and 
can be found in Table [V] An additional uncertainty in the 
scale of the ordinate due to correlated polarization uncertain- 
ties of 2.4% (3.4%) for the points with \x F \ > (x F = 0) is 
not shown. See text for details. 



tion in polarized p + p collisions generated by a gluon 
Sivers TMD would be evidence against large contribu- 
tions from color-octet diagrams for J/ip production. If a 
gluon Sivers TMD is in fact well denned and nonzero, a 
new experimental avenue has been opened up to probe 
the J/ip production mechanism, a long-standing question 
in QCD. As discussed in SectionlH there is a relationship 
between the SSA and the production mechanism of the 
J/ip in the collinear, higher- twist approach, but it is not 
as simple as in the TMD approach. 

Future p + p data from RHIC are expected to improve 
the precision of the current measurement, and a similar 
measurement for J/ip production in SIDIS with a trans- 
versely polarized target could shed further light on the 
production mechanism, as [23| predicts a vanishing asym- 
metry for the color-singlet model in SIDIS but nonzero 
asymmetry for the color-octet model. While no rigor- 
ous quantitative calculations are presently available for 
either collision system, we anticipate that future theoret- 
ical calculations will provide more detailed guidance on 
the implications of the present results. 
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Appendix: Local Polarimetry 

The transverse component of the proton polariza- 
tion is monitored at PHENIX by a measurement of 
the SSA of forward neutron production with a Zero- 
Degree Calorimeter (ZDC) and Shower Maximum Detec- 
tor (SMD), collectively referred to as the local polarime- 
ter. 

The ZDC consists of three modules of hadronic 
calorimeter, 1.7 interaction length each, located approx- 
imately 18 m from the center of the interaction region. 
The SMD consists of two layers of plastic scintillator ar- 
rays located between the first and second modules of 
the ZDC, and it provides horizontal and vertical posi- 
tion information of the detected neutrons. There is an 
additional plastic scintillator in front of the ZDC used to 
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identify and veto charged particles. 

Neutron candidates were selected by triggering on a co- 
incidence between the BBC trigger and a signal in either 
the north or south ZDC. Due to bandwidth restrictions, 
the trigger was prescaled by a factor typically around 100. 
The analysis of the 2006 data used 8.7 x 10 7 events from 
radially polarized collisions, and the analysis of the 2008 
data used 7.3 x 10 7 events from vertically polarized 
collisions. 

Neutrons were selected offline by requiring an energy 
deposit in the ZDC between 20 and 120 GeV along with 
no hit in the veto scintillator (less than 1 minimum- 
ionizing particle). Hits were required in both the hori- 
zontal and vertical planes of the SMD, and neutrons were 
required to be between 0.3 mrad (0.5 cm at the SMD) 
and 1.4 mrad (2.5 cm at the SMD) from the beam. 

The asymmetry is calculated using fits to the az- 
imuthal distribution of the neutrons: 

A(<fr) = An cos(</> — (f>o) (2006, radially polarized) 
A(<fr) = Am sin(</> — <p ) (2008, vertically polarized) 

where 4> is the azimuthal angle relative to vertical, and 



00 is the deviation of the polarization from the nominal 
direction. 

The position resolution of the local polarimeter is ap- 
proximately 1 cm. To determine a systematic uncer- 
tainty on the polarization direction, fits are performed 
for ±1 cm around the center of the azimuthal neutron 
distribution (defined as the position where An is maxi- 
mal). The deviation of the polarization from nominal, in 
radians, was found to be 

0.064 ± 0.040 stat ± 0.086 syst (clockwise) 

0.109 ± 0.038 stat ± 0.036 syst (counterclockwise) 

for the 2006 data and 

0.263 ± 0.030 stat ± 0.090 syst (clockwise) 
0.019 ± 0.048 stat ± 0.103 syst (couterclockwise) 

for the 2008 data. The polarization directions of the 
counterclockwise-going beam in 2006 and clockwise-going 
beam in 2008 were considered to be significant and used 
in the calculation of An for J/tp. 
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